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We show that arrival times for electromagnetic pulses measured through the rate of absorption in 
an ideal impedance matched detector are equivalent to the arrival times using the average flow 
of optical energy as proposed by Peatross et al. [ Phys. Rev. Lett. 84, 2370 (2000)]. We 
then investigate the transport of optical pulses through dispersive media with negative dielectric 
permittivity and negative refractive index choosing the geometry such that no resonant effects come 
into play. For evanescent waves, the definitions of the group delay, and the reshaping delay get 
interchanged in comparison to propagating waves. The total delay times for the evanescent waves 
can be negative in an infinite plasma medium even for broadband pulses. The total delay time is, 
however, positive for broadband pulses in the presence of an interface when the radiation is detected 
outside the plasma. We find evidence of the Hartman effect for pulses when the distance traversed 
in the plasma is much smaller than the free space pulse length. We also show that for a negative 
refractive index medium (NRM) with e(u) = /i(w) the reshaping delay for propagating waves is 
identically zero. The total delay time in NRM is otherwise dominated by the reshaping delay time, 
and for broadband pulses in NRM the total delay time is subluminal. 

PACS numbers: 



I. INTRODUCTION 

Superluminal propagation of pulses has always been 
a fascinating issue for the physics community. Since 
the original analysis of Brillouin and Sommerfeld [l| 
who showed that information transport represented by 
a sharp edge in the wave field can never be superlumi- 
nal, group velocities greater than the speed of light in 
vacuum are often termed abnormal or even meaningless. 
Likewise the time of traverse through a medium or the 
seemingly related arrival time for a wave at a point has 
always been a contentious issue (see [US LI for reviews). 
As waves or wave packets are extremely deformable ob- 
jects, there can be no well-defined start or finish lines for 
them. This has resulted in a variety of time scales that 
have been defined depending on what is being measured 
(see El, IH, IK for some of the approaches). 

The Wigner (group) delay time @ defined as 



dcj)/duj, 



(1) 



where <f> is the phase of the wave, and is based on track- 
ing a fiducial point on a wave packet that moves with the 
group velocity v g = duj/dk, and has been a popular mea- 
sure for the time delay. The Wigner delay time can often 
unsurprisingly be superluminal or even negative as there 
is no causal relationship between the peaks of the input 
and output wave packets J3j- This has often been used 
to term the group delay time as physically unimportant. 
But the group delay time, superluminal or otherwise, has 
been shown to well describe the arrival of electromag- 
netic pulses across absorptive media ; 8j and amplifying 
media particularly for narrow bandwidths and short 
propagation distances. In such cases the apparent su- 
perluminal propagation of the pulse can be explained by 
preferential attenuation or amplification of the trailing 
or leading edge, respectively. However, the experiments 



of Wang et al. [ijl on superluminal pulse propagation in 
an almost nonabsorptive or nonamplifying medium, but 
with highly and anomolously dispersive refractive index 
has been a significant achievement in underlining the im- 
portance of the group velocity in pulse propagation. Now 
ultraslow or ultrafast group velocities including negative 
group velocity have been demonstrated over large dis- 
tances in optical fibers as well Theoretically one 
notes that any pulse with a shape described by a holo- 
morphic function of time can exhibit superluminal prop- 
agation without violating causality or special relativity. 
This is because there is no information in the peak that 
is not contained in the leading edge of the pulse. The 
peak can, in fact, be obtained by a Taylor series expan- 
sion about any point in the leading edge. Information 
can only be encoded by a meromorphic function through 
singularities in the function itself or in some derivative of 
the function. 

For pulses and particularly broadband pulses, Peatross 
et al. ,F2| showed that the arrival time of a pulse at a 
point r can be well described by a time average over the 
component of the Poynting vector S normal to a (detec- 
tor) surface at r as 



<*>r = 



(2) 



Here u is taken to be the unit vector along the normal 
to the given surface. The time of traverse between two 
points (rj,r/) is equal to the difference of the arrival 
times at the two points, and was shown analytically to 
consist of two parts: a contribution by the spectrally 
weighted average group delay at the final point r f 
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u • S(r/, w) [(<9 Re k/du) • Ar] duo 
u -/^°oo S ( r /' w ) dw 
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and a contribution that could be ascribed to the reshap- 
ing of the pulse 

At R = T [exp(-Im k • Ar)E(r,, u)] - T \E(r it u)] , (4) 

which is calculated with the spectrum at the initial point 
IV Here the operator T is 



T[E(r,a,)] 



f°° n e x H*(r,w) 



S(r, u)duj 

(5) 

which represents the arrival time of a pulse at a 
point r in terms of the spectral fields and S(r, u>) = 
Rc[E(r,ix)) x H*(r,w)], represents the Poynting vector. 
Here we take the real parts of the quadratic terms 
since we use complex representation for the fields, i.e., 
e j(k-r-wt) f or a pi ane wave. The total delay time was 
shown to remain subluminal for broadband pulses for 
traversal across a medium with Lorentz dispersion for 
the dielectric permittivity. The most significant aspect 
of this proposal is that it does not involve any pertur- 
bative expansion of the wave number around the carrier 
frequency 

In any experiment, it is not the flow of radiative en- 
ergy that is measured but rather the energy absorbed 
by a detector. Although a proportionality is expected 
between the two quantities, the equivalence is not clear. 
Here we first present our investigations on this aspect 
and conclude that the two delay times are equal for an 
ideal detector with perfect impedance matching. 

We then investigate pulse traversal through dispersive 
media with negative material parameters (e and /i). Neg- 
ative refractive index materials (NRM), which simulata- 
neously have e < and fi < at a given frequency, have 
become exceedingly popular in recent years (see Q 
for recent reviews), particularly since their experimen- 
tal demonstration in artificially structured metamaterials 
EE EE llllll- The phase vector in an isotropic NRM is 
opposite to the Poynting vector and pulse propagation in 
such media can be extremely interesting. One can have 
all combinations of positive or negative phase and group 
velocities in such media |l9j . Pulse propagation in NRM 
has been investigated with a primary focus on the neg- 
ative refraction at interfaces [2(J, |2l], [23 and nonlinear 
NR M I2I I2 I I2H. In the linear regime, Dutta Gupta et 
al. [2|| have studied the group delay time in a cavity 
filled with a NRM in the zero dissipation limit and have 
also shown superluminal propagation in the spectral re- 
gions where the medium only has negative e. Negative 
group velocities in NRM have also been reported [27], [2^ . 

Here we study the delay time for pulse propagation in 
NRM using the time averaged energy flow. We empha- 
size the transport of pulses composed of evanescent waves 
through a plasmalikc medium with e < 0. The traver- 
sal time for evanescent waves itself is a very complicated 
matter j^E Given that NRM can support a host 

of surface plasmon states which can resonantly interact 
with evanescent waves, it becomes even more interesting 



and imperative to study the transport of these pulses of 
evanescent waves. The paper is organized into the fol- 
lowing sections: In Sec. II, we discuss the arrival times 
as measured through the rate of absorption by an ideal 
detector and its connection to the energy flow. The delay 
times for evanescent pulses and pulse transport through 
a plasma (infinite and semi-infinite with a boundary) are 
presented in Sec. III. The delay times for pulse propaga- 
tion in NRM are discussed in sec. IV and we conclude in 
Sec. V with a discussion of our results and their impli- 
cations. 



II. ARRIVAL TIMES AT A DETECTOR 

In an experiment, it is usually the energy absorbed 
by a detector that is observed and not the energy flow 
directly. Here we will investigate arrival time of a pulse 
at a detector by using a time average over the rate of 
absorbed energy in a detector defined as 
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dA(r,t) 
dt 
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dA(r,t) 
dt 



dt 



(0) 



where (dA/dt) is the rate of absorption of energy per unit 
volume inside the detector placed at r. We will call this 
the detector arrival time. The rate of absorption locally 
inside an absorbing medium is given bv [3lll32l| : 



dA 
~dt 



dujdu'[e Lu lm[e(ui)]E*(uj')E(uj) 



This is integrated over the detector volume to obtain the 
total rate of absorption. Obviously the spatial extent 
of the detector is assumed to be small compared to the 
length scales of propagation or spatial pulse widths. Here 
we will assume the pulse propagation along the z axis 
and the detector to be a very thin slab of an absorbing 
medium placed in the path of the pulse. Further the rela- 
tive material parameters (e and [i) are assumed constant 
over the frequency range of interest, i.e., the detector is 
dispersive only over much larger frequency ranges. Here 
we will assume that the absorbing slab has e = (9 + i5) 
and /it = (9 + i5) when the detector slab is impedance 
matched to vacuum, and s = (9 + i5) and fj, = 1 other- 
wise. 

For our calculations we will use a dielectric medium 
with a single resonance Lorentz dispersion: 



e(w) = 



(8) 



and a nondispersive n(ui) = 1. We choose foJi = IOO7 2 
and wo = IOO7. Now consider a pulse of light initially at 
z = whose electric field in time is given by 



E(z = 0, t) = xF. exp[ exp(— iut) 



(9) 
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the arrival times using the average of the Poynting vector. 
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III. ARRIVAL TIMES FOR EVANESCENT 
WAVES 

Now we will consider the arrival times for pulses com- 
posed entirely of evanescent waves. This is analogous to 
the quantum mechanical tunneling of a particle under a 
barrier. Such situations arise directly in the transport of 
radiation across a metal slab or under conditions of total 
internal reflection. Since the phase vectors for the evanes- 
cent waves are imaginary, interesting questions arise re- 
garding their traversal times 0, Q- One of the most 
paradoxical aspects is the saturation of the Wigner de- 
lay time with the barrier thickness - also known as the 
Hartman effect |33j| . 

Now consider the complex wave vector in a medium, 



FIG. 1: (a) The total delay times as a function of frequency 
when measured by using the rate of absorbed energy in a de- 
tector and a time expectation integral over the Poynting vec- 
tor for narrowband pulses, where oj c is the carrier frequency, 
(b) The difference of the arrival time by the detector and the 
arrival time of the pulse using the energy flow equation for an 
impedance matched and a mismatched detector for a narrow- 
band pulse, (c) and (d), respectively, repeat (a) and (b) for 
broadband pulses. 



and the detector slab is at z — c/(107). The pulse 
is propagated forward in time using Fourier transform 
methods and the fields in the absorbing slab are used to 
compute the rate of absorption at subsequent times. The 
effects of multiple scattering within the detector slab are 
included in our calculations. In Fig. ^a), the detec- 
tor arrival times are shown for a narrowband pulse with 
t = IO/7. These parameters are the same as those used 
in Fig. 2(a) of Ref. 0. In Fig-Hfb), we plot the dif- 
ference of this detector time and the arrival time using 
the energy flow given by Eq. (2), and find that there are 
small differences between them particularly near the res- 
onance frequency. However, for an impedance matched 
detector when /j,(oj) = e(w), these differences tend to zero 
as shown in Fig. IHb). This can be traced to the finite 
reflectivity of the detector slab. It is well-known that 
the reflected portion of the pulse will interfere with the 
incoming part of the pulse and can change the traversal 
or delay times. For broadband pulses (r = I/7) also, 
we find that the detector arrival times and the arrival 
times calculated from the average energy flow tend to be 
the same when an impedance matched detector is used. 
We show the detector arrival times and the difference 
between the times in Figs, ^c) and^d). 

Thus we conclude that the detector arrival times and 
the times obtained from the average energy flow of Ref. 
|12| are equivalent for an ideal impedance matched detec- 
tor. In the remainder of this paper, we will only calculate 



k 2 



efi- 



(10) 



In the limit of small imaginary parts of e and /1, one can 
write 



k r = Re(k) ~ \/e r [i r — £iiii — 

c 



ki = Im(k) 



2-\/ E r ^r — SifJ'i C 



(11) 
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where the subscripts r and i indicate the real and imagi- 
nary parts of the quantities. Thus for propagating waves, 
the real part of the wave vector depends primarily on e r 
and [i r while the imaginary part is directly proportional 
to and \ii or the dissipation. This however, becomes 
different for evanescent waves [34]]. To make clear the 
discussion for evanescent waves, we will consider an ab- 
sorbing electric plasma with e r < 0, e, > 0, and fi = fj, r . 
Now, 



(13) 
(14) 



Thus the real part of the wave vector depends on the lev- 
els of dissipation in the medium (ei) and the imaginary 
part of the wave-vector which determines the decay of 
the wave depends on |ey|. This implies, in turn, that the 
definitions of the group delay time and the deformation 
delay time given by Eqs. (3) and (4), respectively, get 
interchanged for evanescent waves. This is an important 
difference for the arrival times of evanescent waves from 
that of propagating waves |5( ■ Note that the same behav- 
ior holds for the case of evanescent waves in total internal 
reflection when the parallel component of the wave vector 
becomes important. This behavior can also be analyti- 
cally continued for larger values of the imaginary parts. 
This essential difference arises because the decay length 
for evanescent waves is determined by e r and not e^. 



A. Pulse traversal in an unbounded plasma 

Now we will investigate the arrival times for evanescent 
pulses inside a plasma. To be specific, we will consider 
the relative dielectric permittivity of the plasma to be of 
the causal form 



G)T=100 Ar=c/100y 



Sl=10 Ar=c/100y 



£ p (w) 



uj(uj + ij) ' 



(15) 



where effectively ujq = and / = 1 in Eq. (8) and the 
relative magnetic permeability /i = 1. For convenience of 
comparison, since the pole of the dielectric function is at 
zero frequency, there is a need to redefine what we mean 
by broadband and narrowband phenomena. To compare 
phenomena across different frequencies, we note that the 
product Or, where uj is the carrier frequency and r is the 
pulse duration in Eq. (9) gives a criterion for broadband 
or narrowband pulses. Typically we take ujt = 1000 or 
100 for narrowband pulses and ujt — 10 for broadband 
pulses. 

We consider that both the source of the radiation and 
the detector are embedded inside the plasma and that 
the plasma is unbounded. This is to avoid any effects 
of scattering from the boundaries and study the inherent 
effects of the plasma on the traversal times. The distance 
between the source and the detector is taken to be Ar. 

In Fig. |3 we plot the delay times obtained for pulse 
traversal inside the unbounded plasma with a plasma fre- 
quency uj p — IO7. First of all, we note that the total 
delay time is negative for a large range of frequencies 
and mostly superluminal below the plasma frequency. 
The total delay time is also dominated by the reshap- 
ing delay time at frequency Q less than u> p . Above the 
plasma frequency for propagating waves, this reshaping 
delay smoothly goes over to the group delay which dom- 
inates the total delay. The dominance of the reshap- 
ing delay below the plasma frequency can be understood 
that inside an infinite plasma, the energy flow for the 
evanescent waves at any point in the dissipative medium 
is determined by the absorbed energy in the region of 
space beyond that point up to infinity. Hence the ab- 
sorption determines the energy transport, and hence the 
delay time as well. In the limit of a monochromatic wave, 
one would expect only this reshaping delay. For larger 
bandwidths, the contribution of the group delay becomes 
appreciable at higher frequencies [uj ~ 77 to IO7 in Fig. 
2(b)]. We also note that the qualitative behavior of the 
delay time does not change appreciably with the increase 
in distance between the source and the detector except 
for the difference in scales (Figs. 2(a) and 2(c), 2(b) and 
2(d), respectively). 



B. Pulse traversal through a bounded plasma 
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FIG. 2: (Color online) The total delay time, the reshaping 
delay time, and the group delay time represented by -, A and 
x as a function of the carrier frequency u in an unbounded 
plasma, (a) Delays for narrowband pulses [lot = 100) and 
Ar =zc/1007. (c) Same as (a), for large distance Ar =zc/107- 
Frames (b) and (d) correspond to broadband pulses (ujt = 10) 
and are similar to (a) and (c). The vertical line drawn at 
cD/7 = 10 separates the propagating waves from the evanes- 
cant waves. Note that the group delay and the reshaping 
delay times interchange their roles for the evanescent waves 
in comparison to the propagating waves. The straight line 
going across the graphs denoted as LL, is the light line for 
free space propagation (At = Ar/c). 



radiation is inside the plasma at a distance Ar from the 
planar interface with vacuum. The detector is taken to 
be in vacuum just outside the interface. This is more 
physical because there would be an interface (impedance 
mismatched) involved with the detector anyway. A cor- 
responding physical situation would be an atom located 
within the plasma and emitting radiation, whose leakage 
is detected outside the plasma. For simplicity, we con- 
sider only waves with a zero parallel wave vector. Then 
the radiation does not couple to any surface plasmon 
modes of the plasma-vacuum interface. Thus our source 
would be an infinite shield of current parallel to the inter- 
face. The interface has an important effect of changing 
the amount of energy that reaches the detector via the 
transmittance of the interface given by 



2^7 



(16) 



where T p is the Fresnel transmission coefficient for the 
P-polarized light. Ef — 1 is the relative dielectric per- 
mittivity in free space. Thus the field at the detector is 
given by 



H(r / ,a;)=T p (a;)e lk ' Ar H(r i ,a;), 



(17) 



Here we will consider the situation of a plasma with a 
semi-infinite extent. We will consider the source of the 



where k is the wave vector in the plasma. Note that the 
final delay times however do not depend on the polar- 
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FIG. 3: (Color online) The various delay times plotted for 
the bounded semi-infinite plasma. As in Fig. 2, the graphs 
(a) and (c) correspond to the narrowband pulses and (b), (d) 
correspond to the broadband pulses for different distances. 
The symbols are similar to Fig. 2. 



ization at normal incidence(k|| — 0). The geometry of 
a semi-infinite plasma that we consider here also avoids 
any coupling with the slab resonances such as Fabry- 
Perot resonances or surface plasmon polaritons. Thus, 
we can study purely the effects of the intrinsic plasma on 
the traversal time. 

In Fig. |3 we plot the delay times for pulses from a 
source within a plasma for different bandwidths and dis- 
tances of the source from the boundary. First of all, we 
note that the reshaping delay time is negligible compared 
to the group delay time even for Q < u> p . The group de- 
lay dominates the total delay time for small source to 
interface distances (Ar = zc/107 and ZC/IOO7). This im- 
portant difference from the case of an unbounded plasma 
results because the presence of the boundary causes a 
reflected evanescent wave. Now energy transport is pri- 
marily determined by the phase difference of the incident 
evanescent wave and the reflected wave, and not only by 
the dissipation in the medium. Thus the group delay time 
plays the determining role. Secondly it should be noted 
that the total delay time is almost always positive except 
near the plasma frequency for narrowband pulses and 
small Ar (=zc/107 and =zc/1007). One notes that neg- 
ative delay times result for narrowband pulses (iDr = 100) 
and small Ar even for propagating waves (e > 0) near 
u> = u>p. This is a consequence of the k = mode at 
e = 0. This negativity goes away for larger bandwidths. 

At a much larger source to interface distance 
(Ar =zc/7), the deformation of the pulse contributes 
appreciably to the total delay time. We show the delay 
times in Figs. 4(a) and 4(b) for Ar —zc/j at different 
bandwidths. The behavior of the reshaping delay time 
(lu < ujp) tends to that of the behavior in an infinite 



FIG. 4: (Color online) (a) The various delay times plotted 
for large source-boundary distance (Ar —zc/^y) and narrow- 
band pulses (pT = 100) for a bounded semi-infinite plasma. 
The symbols are similar to Fig. 2. (b) Same as (a), but for 
broadband pulses (cut = 10). The Hartman effect can be ob- 
served in graphs (c) and (d) which correspond to the total 
delay times for various parameters shown. 



plasma while the group delay time strongly moderates 
this contribution to the total delay, and the total de- 
lay time is positive almost everywhere. Surprisingly we 
note that there is a small region of frequencies where 
the total delay time goes negative even for broadband 
pulses, and at this large distance involved. However, we 
note that the spectral width of the region where the to- 
tal time becomes negative, reduces with increasing pulse 
band width. Hence in the limit of very large bandwidths, 
we expect this spectral width to go to zero asymptoti- 
cally. 

We note the presence of a Hartman effect in our calcu- 
lations as well. In Figs. 4(c) and 4(d), we plot the delay 
times with respect to carrier frequency without scaling 
with respect to the temporal pulse width. We find that 
over a large range of carrier frequencies below the plasma 
frequency, the delay time is almost the same for various 
distances involved (Ar =zc/107 and ZC/IOOO7). This is 
seen for both broadband pulses as well as narrowband 
pulses for these Ar. This saturation of the delay time 
with distance is a generalization of the Hartman effect for 
broadband pulses which is usually noted for monochro- 
matic evanescent waves. However, for much larger dis- 
tances (Ar =zc/j) which are comparable to the spatial 
pulse width in free space (I = cr), the deformation takes 
over and the Hartman effect is lost. It is important to 
mention that in our case it is the dissipative nature of 
the medium that destroys the Hartman effect. Finally, 
we note that, the negativity or superluminal delays are 
always a fraction of the pulse width, and thus does not 
imply any violation of causality in all these cases. 



IV. ARRIVAL TIMES IN NEGATIVE 
REFRACTIVE INDEX MEDIUM 

In this section, we will discuss the arrival times for 
pulses propagating through negative refractive index 
medium. A sufficient condition for negative refractive 
index is e < and /i < at any frequency [l4| . 

First, we note that an unusual situation arises when 
the value of e becomes equal to that of /j, (a case of prop- 
agating waves). For this case, the reshaping delay turns 
out to be zero. Note that the value of the refractive index 
becomes equal to e or /i when e — fi. The electric field 
at the initial point is given by 
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In an unbounded medium, the electric field at the final 
point is related to that at the initial point as 
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FIG. 5: (Color online) Dispersion of the real and imaginary 
parts of e and given by Eqs. (15) and (23). The vertical 
lines at frequencies 57, 9.4347, and IO7 separate the frequency 
regions where the field modes are either evanescent or propa- 
gating. See text after Eq. (23). 



where r/ = + Ar and k is the wave vector in the 
medium. The magnetic field is related to the electric 
field through the Maxwell's equation, 



H(r;,w) =y— -= re 
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Using these, we calculate the delay time for pulse prop- 
agation between the initial and the final positions and 
see that the total delay time consists of only one nonzero 
term which is the group delay time, 



At = 
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This means that the reshaping delay time identically 
vanishes and no reshaping of the pulse takes place in 
a medium with s = fx. Even for a bounded medium 
with an interface separating the given medium from vac- 
uum, due to perfect impedance matching, we note that 
the transmission coefficient through the interface is unity. 
This means that the dispersion in the transmission of the 
pulse through the interface plays no role and the reshap- 
ing delay disappears here too. 

For concreteness, we consider the following causal dis- 
persive forms for e and /1: The Drudc Lorentz form for 
e given by Eq. (15) and a Lorentz dispersion for [i given 
by 



J 0m 



(23) 



For convenience, we take u>o m = 57, uj^ n — 647 2 , and 
the rate of dissipation 7 in [i to be the same as that 
in e. This results in an electric plasma [e < 0, fj, > 
0) for < oj < 57, a negative refractive index medium 
(e < 0, fj, < 0) for 57 < u> < 9.4347, an electric plasma 
(e < 0, fi > 0) for 9.4347 < uj < IO7, and a positive 
refractive index medium (e > 0, fi > 0) for lj > IO7. 
The dispersions of e and fj, are shown in Fig. 5. 

First, we study the behavior in an infinite medium. We 
plot the delay times (total, group and reshaping delays) 
for a pulse in Fig. 6. We note that at low frequen- 
cies, the delay time is negative as in an infinite plasma. 
However, there is a large peak in the total delay time at 
u> ~ 37. The group delay almost exclusively contributes 
to this. This is presumably due to the rapid increase in 
Ke(fi). In the negative refractive region, when there are 
propagating waves in the medium, there is appreciable 
contribution from both group and reshaping delays. The 
reshaping delay time is mostly negative. We note that 
the total delay time is always positive as well as sublu- 
minal in the negative index frequency region. 

Now, we take a semi-infinite medium with the above 
material dispersion. As in Sec. Ill B, we take the source 
to be inside the semi-infinite medium at a distance Ar 
from the interface with the vacuum and detect the radi- 
ation outside. As before, we consider only waves with a 
zero parallel wave vector (k|| =0). In Fig. 7, we plot the 
delay times for different distances in both broad and nar- 
rowbands. At low frequencies for narrowband pulses, we 
have the behavior in a plasma [Figs. 7(c) and 7(d)]. Even 
for frequencies, ui > IO7, when the index is positive, the 
nonunit dispersive magnetic permeability (/z < 1) gives 
rise to a large reshaping delay for broadband pulses. We 
note that the reshaping delay plays a major role in deter- 
mining the times in the negative refractive index region. 
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FIG. 6: (Color online) The various delay times plotted as a 
function of to in an unbounded medium which behaves like 
a plasma and with negative or positive refractive indices for 
certain frequency ranges for broadband pulses (lot — 10) and 
Ar =zc/107- The vertical lines at frequencies 57, 9.4347, and 
IO7 separate the frequency regions where the field modes are 
either evanescent or propagating. The symbols are similar to 
Fig. 2. 
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FIG. 7: (Color online) (a) The various delay times plotted 
as a function of Q in a bounded medium. The parame- 
ters and symbols are similar to Fig. 6. (b) Same as (a), 
with Ar --ZC/TOO7. Frames (c) and (d) correspond to nar- 
rowband pulses (lot — 100) with source-boundary distances 
Ar =zc/107 and Ar --ZC/TOO7. 



In fact for broadband radiation, the bulk of the total 
delay time comes from the reshaping delay [Figs. 7(a) 
and 7(b)]. The group delay time disperses rather vio- 
lently near Co = LOQ m . For narrowband pulses {lot = 100), 
the total delay time can become negative near regions 



where a change of sign occurs in n or e. But for broad- 
band pulses (lot = 10), this negativity disappears even 
for short distances of propagation. We also note that 
for broadband pulses, the total delay time can be signif- 
icantly larger (of the order of a pulse width) and it is 
always subluminal when the refractive index is negative. 



V. CONCLUSIONS 

Pulse propagation through dispersive materials has al- 
ways raised counterintuitive issues regarding superlumi- 
nal propagation and negative delay times. Measuring 
the arrival times of electromagnetic pulses via the aver- 
age flow of energy was first proposed by Schwinger [3{| 
and the neat decomposition of this delay time into the 
average group dela y ti me and the reshaping delay time 
by Peatross et al. |T2| has made it an attractive candi- 
date for describing the traversal of pulses. In fact these 
times have also been exp erimentally measured underlin- 
ing their importance |3a |. 

In this paper, we have first shown that the arrival times 
based on the averge energy flow arc equivalent to the 
times measured by the rate of absorption in the detec- 
tor volume for an impedance matched perfect detector. 
Thus the intuitive feeling that the two times should be 
equivalent has been confirmed by rigorous comparison. 
An impedance mismatch, however, renders the two times 
nonequivalent. 

We then investigated the delay times using the average 
energy flow for evanescent pulses and demonstrated an 
important difference from the propagating pulses. The 
very definitions of the average group delay time and the 
reshaping delay time for the evanescent pulses get inter- 
changed due to the (primarily) imaginary wave vector 
for evanescent waves. Thus <9Re(k)/9a-> contributes now 
to the deformation of the pulse via the dissipation and 
Im(k) contributes to the group delay. 

We then evaluated the delay times for pulses traversing 
a plasma with negative s(lo). We have shown that in an 
infinitely extended plasma, the delay time is primarily 
determined by the reshaping delay time and is usually 
negative. This is because the energy flow at any point in 
the infinite plasma is essentially determined by the dis- 
sipation in the regions up to the point of detection. On 
the other hand, in a bounded plasma when the radiation 
is detected outside the boundary (in vacuum), the group 
delay time dominates and the total delay times are usu- 
ally positive and subluminal for large enough frequency 
bandwidths associated with the pulses. The interface 
also modulates the delay times via the dispersion of the 
transmission coefficient. The reflected evanescent waves 
play an important role in the energy transport. We also 
note the Hartman effect in the context of energy trans- 
port for evanescent pulses when the source to boundary 
distance is small compared to the free space pulse length. 

In the case of negative refractive index materials, mod- 
eled by causal e(lo) and fJ-(i-o), the total delay times are 
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dominated by the reshaping delay times. The group de- 
lay time contributes largely near the transition frequen- 
cies when Re(e) and Im(/z) change signs. We have proven 
an important result that the reshaping delay time is iden- 
tically zero for a medium with e(w) = (J>(u>). 

We should point out that the delay times that we have 
calculated here are due to the intrinsic dispersive nature 
of the material parameters, e{ui) and The geome- 

tries that we have chosen, that of infinite or semi-infinite 
media containing the source ensure that no geometri- 
cal resonances such as Fabry-Perot resonances for slabs, 
are present to affect the times calculated. Although the 
media with negative e(w) and /Lt(w) can support a vari- 
ety of surface plasmon resonances even on the surface of 
semi-infinite media, our calcuations for normally incident 
waves with k|| = ensures that radiation cannot couple 
to these resonances. 

Finally we would like to point out that although we ob- 



tain superluminal or negative total delay times in many 
cases, for example in a plasma, the superluminality or 
the negativity is always a fraction or of the order of the 
pulse width. Thus for a Gaussian pulse with an infinite 
support (in principle), there is no violation of causality 
implied by our results. This is in spite of our calculations 
based on energy transport. 
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